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Abstract 

The bottomonium spectral functions at finite temperature are analyzed by employing QCD sum rules with the 
maximum entropy method. This approach enables us to extract the spectral functions without any phenomenological 
parametrization, and thus to visualize deformation of the spectral functions due to temperature effects estimated from 
quenched lattice QCD data. As a result, it is found that T and ijt survive in hot matter of temperature up to at least 
2.3T C and 2.1T C , respectively, while Xw and^/,1 will disappear at T < 2.5T C . Furthermore, a detailed analysis of the 
vector channel shows that the spectral function in the region of the lowest peak at T — contains contributions from 
the excited states, T(25) and T(35), as well as the ground states Y(1S). Our results at finite T are consistent with the 
picture that the excited states of bottomonia dissociate at lower temperature than that of the ground state. Assuming 
this picture, we find that T(25) and T(35) disappear at T = 1.5 - 2.0T C . 
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1. Introduction 

Properties of high density and temperature matter is one of the most exciting subjects of hadron physics. Quantum 
chromodynamics (QCD) predicts matter composed of quarks and gluons in the form of unconfined plasma phase, 
called quark gluon plasma (QGP). In such a phase of matter, production of quarkonium, a bound state of a heavy quark 
and a heavy antiquark, is expected to be suppressed. Especially, J/tfr suppression has traditionally been considered as 
a signal for the formation of QGP JH01 m high-energy heavy ion collisions. There, the charmonia including J/ift may 
dissociate in the QGP due to temperature effects such as the color Debye screening. Such suppressions were observed 
in the heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC). A similar phenomenon may occur in the 
bottomonium spectrum. Recent report [3j has suggested significant modification of the bottomonium spectra from 
the comparison between the P-P and Pb-Pb collisions at the Large Hadron Collider (LHC). The data indicate larger 
suppression of the excited states than the ground state. One purpose of the present paper is to study this phenomenon 
by using QCD sum rules with MEM. It is thus interesting to see whether there are qualitative/quantitative similarities 
and differences between the behavior of the charmonia and bottomonia spectra at finite temperature. 

In order to quantify and predict when and how the quarkonium spectrum is modified, various theoretical ap 



proaches have been developed. They include Lattice QCD [0, El S 0. 



eloped, iney : 
iKHEl, AdS/QCD 01j2(jj2l]j2^ 



BE 



121 [LJ, QCD sum rulesj 14U15J 
23||24|, and effective field theories II251 12611 as 
Most of the previous studies were devoted to the 



well as potential modelslTiriHlElii [H El HlHlf 
charmonium spectra, while, for bottomonia, there are only a few theoretical results 
and it is important that they are examined by independent methods 
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In a previous work 111 811 . we have investigated the behavior of charmonium spectral functions at finite temperature 
from QCD sum rules with the maximum entropy method (MEM). Our approach enables us to extract directly the 
shape of the spectral functions. Therefore, this is a suitable tool to study the deformation of the spectral functions 
upon the change of the temperature. In this paper, we will study the bottomonium spectrum by using the same method. 
We point out that the bottomonium spectral function is crucially different from that of charmonium in the sense that 
there are several excited states below the continuum threshold. For example, T has the ground state (IS) and the 
excited states (25, 3S). 

This paper is organized as follows: In Section [2] we discuss QCD sum rules for heavy quarkonia at finite tem- 
perature and demonstrate our method to extract spectral functions with MEM. In Section [3j we show the obtained 
spectral functions at zero and finite temperature. Here, we also investigate the behavior of the excited states of the 
bottomonium vector channel. Section|4]is devoted to the summary and conclusion. 



2. Analysis procedure 

This section overviews our analysis method of QCD sum rules with MEM. Historically, QCD sum rules for heavy 
quarkonia were introduced in ll38l l39ll and then elaborated in lf40h . Their Borel transformation were calculated in [ 41 1 . 



Applying MEM to QCD sum rules was successfully performed for the p meson in 114211 and the nucleon in (43 



We follow the same procedure as 01811 for the QCD sum rules of heavy quarkonia at finite temperature. 



2.1. QCD sum rules for bottomonia 

Let us start with the current correlation function 

n J (q) = ij ( / , w- , " > <y|./'(v) / /: (0)|>. (i) 

where J stands for the pseudoscalar (P), vector(V), scalar (S), and axial-vector (A) channel. Each current is defined 
as j p = bjsb, ]Y = by M b, j s = bb, and jfi = (q^qv/q 2 - g^jby^b with b being the bottom quark operator. For the 
axial-vector current, the projection to the transverse components for the^i states is carried out. At finite temperature, 
the currents of the V and the A channels generally have two independent components. We assume spatial momentum 
of the bottomonia to be zero so that only one component becomes independent : qf 1 = (a), 0). Then, one can define the 
dimensionless correlation functions as f[ P S (q 2 ) - H p,s (q)/q 2 and tl v ' A (q 2 ) = Iff, ' (<?)/(— 3q 2 ). 

Using the operator product expansion (OPE), one can expand the operator j J (x) j J ~(Q) of Eq.([TJ as a series of local 
operators O n with mass dimension n. Then, the dimensionless correlation functions are given as 

nV) = 2 c '^ 2)<0 " ) - (2) 

n 

If the scale of the gluon condensates is smaller than the separation scale: Am 2 b - q 2 » (G) ~ (Aqcd + aT + bu) 2 , one 



can assume that all the temperature effects are included in the expectation values of the local operators (O n ) 1451, 1 1 711 . 
Thus the Wilson coefficients C/ (q 2 ) can be considered to be independent of T, as long as the temperature is not too 
high. 

In order to improve the OPE convergence and suppress contributions of higher energy states, we perform the Borel 
transformation on Eq.(f2). Then, the correlation function can be written down as 

M J (v) = ne- v A J (v)[l + a s (v)a J (v) + b J (v)<f> b (T) + c J {v)4> c {T)], (3) 

where we use the dimensionless parameter v = Am 2 IM 2 with M standing for the usual Borel mass. Now the condition 
for the separation scale is given by d\M 2d » {G d ) 11711 . The first term is the leading order of the OPE corresponding 
to the free current correlation function. The second term stands for the perturbative a s correction. The third and 
fourth terms include contributions of the scalar and twist-2 gluon condensates of mass dimension 4. Here, <ph and <p c 
are defined as 

** (r) = o7?W Go(:r) ' (4) 

9(4f7z, Y 



<Pc(T) = T7r ^G 2 (r), (5) 



4tt 2 
3(4m2)2 

where Go(T) and G2(T) are the scalar and twist-2 gluon condensates at finite temperature. Go and Gi are defined as 
G (T) = <f G^G"^} and (w^V - i^ v )G 2 (r) = (f G a /G avp ), where ^ is the four velocity of the medium. In our 
previous work on charmonium, we have added the scalar gluon condensate of mass dimension 6 which was found to 
be small. Thus, we can safely neglect it because its coefficient is strongly suppressed by the bottom quark mass. The 



detailed expressions of the Wilson coefficients A J (v), a J (v), b J {v) and c J (v) are given in Ref. IU7I1 



The temperature dependences of the gluon condensates are obtained from the approach proposed in Refs. I1141I15I1 . 



where the dimension-4 gluon condensates are related to the energy-momentum tensor, which can be expressed in terms 

'0 



of the energy density e, the pressure p and the strong coupling constant a s . In concrete, Gq(T) = GT e - f^[e(T) 



3p(T)] and G 2 (T) = -^[efT) + p(T)]. We then utilize the results of quenched lattice QCD |46J,|47|] to obtain the 



temperature dependence of p(T) and a s (T). 

It should be noted here that, in the quenched approximation, the value of the critical temperature T c is about 260 
MeV l46Tl . while the cross-over temperature of full QCD is estimated to be in the region of 145-165 MeV [48, 49 ]]• 



In order to make the predictions more realistic, full lattice QCD with physical quark masses 115 111 may be applied. 
We, however, note that the gluon condensates are not directly extracted from energy density and pressure in full QCD 
because of the light quark contribution. Therefore, in the present approach, we have to assume that Go{T), G2(T) 
and a s (T) are functions of T/T c . For the scalar gluon condensate, this approximation is known to be a good around 



T/T c ~ 1.0 152111511 . Thus this approach may provide a good approximation for charmonium systems. As we shall 



see, however, bottomonium systems receive significant modification at higher temperatures where not only Go(T), 



unknown for full QCD, becomes more important [ 17) but also Gq(T) in full QCD deviates from the scaling behavior. 
Therefore, resultant temperature dependence of the spectral functions, which will be discussed below, provides no 
more than qualitative guides for the full QCD case. Nevertheless, we emphasize that the results obtained in this 
paper, rather than being entirely applicable to full QCD, can provide useful information on the difference between the 
behavior of the S-wave and P-wave (or charmonium and bottomonium) states. 

The sum rule is constructed from the dispersion relation derived from the analytic properties of the correlation 
function of Eq.([TJ. For the spectral function p J (a>), we obtain 

Jr»oo 
dco 2 e- m2/4 <p J (aj), (6) 
o 

where the left-hand sides are equal to Eq.(|3). Note that, in the vector, scalar, and axial-vector channels, there is an 
additional constant term to Eq.Q at finite temperature, which originates from a pole at a> — in p J (oj) and is called 
scattering term [53]. Although the contribution of this term, proportional to e~ m <>l T , should be much smaller than the 



charmonium case, we differentiate Eq.© with respect to v to eliminate this contribution: 

M J (y) = - da> 2 aj 2 e-™ 2/4m lp J (<o). (7) 

av 4m- b Jo 

To compare all the channels on the same basis, we also differentiate the sum rule of the pseudo-scalar channel. The 
validity of this procedure in the heavy quark sum rules was discussed in Ref. 15411 . Also, we have checked that the 



same results are obtained from both the original and differentiated sum rules for the pseudo-scalar channel. 

2.2. MEM analysis of QCD sum rules 

In this subsection, we demonstrate our method to extract the spectral function p J (a>) fromEq.(|7]i. The conventional 
methods of analyzing QCD sum rules assume a particular form for the spectral function, the most popular one being 
the "pole + continuum" form. By contrast, such an assumption is not necessary in our method [42] as the shape of 
spectral functions is directly extracted from the MEM. 



'Whereas the phase transition in pure SU(3) theory is a first order one, full QCD at physical quark masses exhibits a cross-over jsoh - 
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Let us now briefly summarize the procedure of the MEM analysis. The basic idea is Bayes' theorem : 



F P[M\H] 

where p and M are corresponding to the spectral function and the OPE in Eq. (O, respectively. H denotes prior 
knowledge on p such as positivity and its asymptotic values. P[p\MH] represents the conditional probability of p 
given M and H. On the right-hand side, P[Ai\pH] is called the "likelihood function", and P[p\H] stands for the "prior 
probability". P[A1|fl] is only a normalization constant and can be ignored in later discussion since it does not depend 
on p. In order to maximize P[p\AiH], we estimate P[M\pH] and P[p\H]. 
The likelihood function is written as 

P[M\pH] = e~ m , (9) 
rn 1 p» lM(v)-M p (v)] 2 

Up] = 7T, 7 dv j— . (10) 

2(v max - v min ) J Vmjn cr 2 (v) 
Here, M{v) is obtained from the results of the OPE and corresponds to the left-hand side in Eq.©, while M p {v) is 



defined as the right-hand one. cr(y) stands for the uncertainty of M(v) (see Ref. 14211 ). On the other hand, the prior 
probability is written as 

P[p\H] = e aSlp] , (11) 
p(u>)' 



-r 

Jo 



S[p]= I dcu) 



p{(jS) - m(u>) - p(u>) log 



(12) 



\ m{u>) j 

where S [p] is known as the Shannon-Jaynes entropy and a is introduced as a scaling factor. m(u>) is called the default 
model and determines the spectral function when there is no information from the OPE. For the default model, we use 
a constant corresponding to the perturbative value of the spectral functions at high energy. 
Using Eq.© and Eqs.dTTT). one can rewrite Eq.® as 

P[p\MH] oc P[M\pH]P[p\H] 

= e Q{p \ (13) 
Q[p] = aS[p]-L[p]. (14) 

In order to determine the most probable p(u>), we find the maximum of the functional Q[p] by the Bryan algorithm 



I 55H . It can be proven that the maximum of Q[p] is unique if it exists, as it is shown explicitly in \5(. 

Before starting the MEM analysis, we have to determine the criterion for the range of the dimensionless parameter 
v which is used for the investigation. We determine the upper bound v max from the criterion that the perturbative a s 



correction term contributes less than 30 % of the leading term as in H17H . The reason why we use the perturbative part 
rather than the term of the highest gluonic condensate as a criterion of the OPE convergence, is that for bottomonium 
the terms proportional to the gluon condensates are strongly suppressed due to the large bottom quark mass, as seen 
in Eqs.fHJi and (Q. Therefore, there is a region of v, where the a s correction is forbiddingly large although the 
contributions of the condensates seem to converge. On the other hand, for the lower bound v m j n , we do not have to 
impose such a criterion unlike in usual QCD sum rule analyses, since we are not assuming the pole dominance of the 
dispersion integral. Thus, we choose the value which obtains the highest peak as the criterion of v m ; n . The range of v 
(namely, Borel window) for each channel is summarized in Table Q] We have checked that the obtained results do not 
depend on the chosen values of these bounds. 



Channel 


Vector (T) 


Pseudoscalar (?#,) 


Scalar (xbo) 


Axial-vector (xb\) 


^min ^max 


4.00-8.23 


4.30-9.33 


4.00-7.33 


3.50-5.83 



Table 1 : Range of dimensionless parameter v. 
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3. Results and discussion 



3.1. Analysis of mock data 

In order to estimate the resolution of the MEM, we should carry out a test analysis by using mock data generated 
from the experimental information for the vector channel, delta functions for the bound states and a smooth step-like 
function for the continuum: 

Imn v (i) = V — • T s (e + O— 7 • 6(s - ml) + continuum, (15) 

res 

where, the sum is taken for the these bottomonium states, T(15), T(25) and T(3S). tor represents their masses, 
rg(e + e~) is the corresponding decay width to an electron-positron pair, and a is the fine structure constant (see Ref. 
J40D). We construct mock data by substituting the phenomenological spectral function into Eq.O and evaluating the 
integral numerically over s. We also consider the mock data for the individual bottomonium states, where a single 
state + continuum is taken for Eq.dTSb. 

The results of the mock data for the individual bottomonium states are shown in Fig. [T]as the dotted lines. The 
result for the full spectrum with all the bound states is shown in Fig. Q~|as the solid line. One sees that the individual 
peaks are not resolved in the full spectrum. Namely, they are combined into a single peak after the MEM analysis. 
Moreover, the residue of the single peaks is consistent with the sum of the residues of T(15), T(2S) and Y(3S). It 
should be noted that the peak position of the combined spectrum is located at a higher energy than the T(15) peak. 
The values of the peak positions and the residues are summarized in Table [2] Later, we will use this property of the 
combined peak, when we extract the temperature dependence of the excited bottomonium states. 

0.07 
0.06 
0.05 

_ 0.04 

3. 

°~ 0.03 
0.02 
0.01 


6 8 10 12 14 16 18 20 
to[GeV] 

Figure 1: Spectral functions extracted from mock data based on Eq. )15t with MEM. 





T(15) 


T(25) 


T(35) 


Y(1S +25 + 3S) 


mass (exp.) [GeV] 


9.46 


10.02 


10.36 




residue (exp.) [GeV] 


0.0270 


0.0123 


0.0089 


0.0483 


mass (from mock data) [GeV] 


9.46 


10.06 


10.42 


9.61 


residue (from mock data) [GeV] 


0.0303 


0.0170 


0.0147 


0.0506 



Table 2: Mass and residue values obtained from mock data based on Eq. )15) with MEM. 



mockdata(1S+2S+3S) 




« mock data(1S) 




1 mock data(2S) 




mockdata(3S) 


I 
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3.2. Bottomonia at zero temperature 

The results of the MEM for the spectral functions of the vector, pseudoscalar, scalar and axial-vector bb systems at 
zero temperature are shown in Fig. |2] We employ the bottom quark mass fhb(mi,) = 4.167 + 0.01 3GeV [57], the strong 
coupling constant a s (M z ) = 0.1184 + 0.0007 Q and the vacuum gluon condensate G™ c = 0.012 + 0.0036GeV 4 
lT38h . Each spectral function shows a clear peak at around to ~ lOGeV. The average peak height and the estimated 
error of the average height in the MEM procedure are given by the horizontal and vertical lines. The lines are 
drawn at the range of to where the average is taken. One sees that the lowest-energy peaks of the S-wave channels 
are statistically significant. On the other hand, the lowest P-wave peaks have heights of the same order as their 
error-bars, the scalar channel being a bit higher while the axial-vector channel is lower. This means that with the 
precision presently available, we cannot make any strong statement about the existence of these peaks and about 
their behavior at finite temperature. Therefore, the results about the lowest P-wave peaks obtained in this paper 
should not be considered to be fully conclusive. Furthermore, for all peaks appearing at higher-energy, their statistical 
significances are not good. The lowest peaks are located at mr = 9.56GeV, m nb = 9.51GeV, m Xm = 10.15GeV 
and m Xbi = 10.42GeV, respectively. These values are somewhat higher than the experimentally observed masses 
(w exp . = 9.460GeV, 9.389GeV, 9.859GeV, and 9.893GeV, respectively). In fact, they are consistent with our analysis 
of mock data since the obtained peaks contain contributions from both the ground and excited states and their positions 
are shifted to higher energies. Validity of this picture is confirmed by evaluating the residue of the peak for the vector 
channel and comparing it with the residue obtained from the leptonic decay width (the value obtained from the present 
analysis is 0.0476 GeV, which should be compared to Table|2]i. 




Figure 2: Spectral functions of bottomonia at zero temperature. Upper left : vector (T), upper right : pseudoscalar (rj),), lower left : scalar (xim), 
lower right : axial- vector (xi,i ). Vertical lines stand for statistical error of peaks obtained from MEM, while horizontal lines represent corresponding 
range of ui, for which the error is calculated. 
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We have checked that such contributions of the excited states are also present in analyses by usual QCD sum rules, 
where the "pole + continuum" assumption is used, i.e. ground state cannot be separated as a single pole. 

3.3. Bottomonia at finite temperature 

The results of the spectral functions at finite temperatures are shown in Fig. [3] All the channels show the same 
qualitative behavior. First, the peak position undergoes a shift to a lower energy with increasing temperature. Next, 
the peak gradually lowers, becomes broader and simultaneously shifts to slightly higher energies. At the same time, a 
continuum-like structure grows in the low energy region, penetrates into the peak regions and moves downward. As 
a further point, let us mention the crucial difference between the behavior of S-wave and P-wave channels. Quantita- 
tively, the vector and pseudoscalar bb states, T and T]b, remain as clear peaks up to T/T c = 2.0 and may still survive at 
T/T c = 2.5. In the case of T, one sees a bump even at T/T c = 3.0. On the other hand, XbO ® a <&Xb\ seem to disappear 
at T/Tc = 2.0 - 2.5. 

Since the uncertainties involved in the calculation are quite large, as indicated by the error bars in Fig. |4] we 
can presently not make statements about specific numbers for the melting temperatures of the investigated states, but 
can only restrict the ranges of temperatures within which the peaks disappear. Concretely, we define the range of 
the melting temperatures as follows. The upper limit is determined as the temperature where the bump (extremum) 
disappears, while the lower limit taken as the temperature where the error bar exceeds the lowest-energy peak height of 
the spectral function, namely when the peak ceases to be statistically significant. The error bars for each temperature 
are shown in Fig. [4] The resultant upper and lower limits of the dissociation temperatures are summarized in Table 
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Figure 3: Spectral functions of bottomonia at finite temperature. Upper left : vector (T), upper right : pseudoscalar lower left : scalar ixm), 
lower right : axial-vector (Xbl)- 
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Channel 


Vector (T) 


Pseudoscalar (rji,) 


Scalar (xbo) 


Axial-vector {xb\ ) 


T/T c 


>2.3 


> 2.1 


1.3-2.5 


< 2.5 



Table 3: Upper and lower limits of the dissociation temperature ranges for the lowest bottomonium states. The precise definition of these limits is 
given in the text. 



[3] If we were able to calculate the OPE data with better precision, the lower (upper) limit of the temperature ranges 
would be increased (decreased). 

In order to confirm that the obtained results are caused by genuine physical effects, we have to check possible 
contributions of MEM artifacts at finite temperature. First, as the contributions of the gluon condensates increase 
at finite temperature, their uncertainties magnify the OPE error. Therefore, it is expected that the resolution of the 
MEM is reduced and the peaks of the extracted spectral functions become broader. Thus, to investigate this effect, 
we reanalyze the spectral functions by using the OPE at T — with error at finite temperature. The result for each 
channel is shown as dashed lines in Fig. [4] As one can see, although the heights of the peaks are reduced partly due 
to the MEM artifact described above, the peaks are still present in the S-wave (P-wave) channels with the error of 
T/T c = 3.0 (T/T c = 2.5). We also stress that the MEM artifact does not shift the peak position. From this analysis, 
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Figure 4: The solid lines and their error bars show the spectral functions of bottomonia obtained from finite temperature OPE data. The dashed lines 
stand for spectral functions at zero temperature using error at corresponding finite temperature. Upper left : vector (T), upper right : pseudoscalar 
{rib), lower left : scalar (Xbo), lower right : axial- vector (Xbl)- 
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we conclude that the disappearance of the peaks at the finite temperatures is caused by physical effects and not due to 
an MEM artifact. 



3.4. Excited states of bottomonia 

In Section 3.1, we showed that the spectral function extracted by the MEM contains contributions of the excited 
states in the lowest peak. In order to extract finite temperature effects on the excited states from the spectral functions, 
we analyze the residue of the lowest peak of the vector channel. However, one cannot naively integrate the spectral 
function in the region of the peak because the spectral function is contaminated by the continuum, which is not 
negligible particularly at high temperatures. 

In order to exclude the continuum contributions and to estimate the sum of the residues only of the ground and 
excited bb states, we fit the obtained spectral functions by using a Breit-Wigner (or Gaussian) function for the peak, 
and the continuum parametrized in the form of the leading order perturbative result. Specifically, we take 



/(*) = 



\A\ 2 



r 



2n (x - m) 2 + T 2 /4 



1 

8^? 



1 - 



4a 2 



2 + 



4a 2 



(16) 



for fitting the MEM results. The four fitting parameters, \A\ 2 , m, F, and a correspond to the residue, peak position, 
width, and continuum threshold, respectively. Note that a coincides with m\, in the perturbative calculation. These 
parameters are fitted by the Levenberg-Marquardt method ll59l l60ll . Furthermore, in order to exclude a possible initial 
value dependence of the fitting procedure and to investigate the existence of local minima, we take 200 initial values 
generated randomly for the four fitting parameters at each temperature. 

The obtained residues of the T peak, |/i| 2 , with increasing T/T c are plotted in Fig. [5] The left panel shows the 
results of the fitting with the function Eq.([T6ll. For each temperature, 200 results corresponding to initial values are 
plotted. For some temperature, we find multiple solutions, which are supposed to be local minima solutions of the 
least-square function^ 2 in the L-M method. In the case of the Breit-Wigner + continuum fitting, we find that the local 
minimum form three clusters, top, middle, and bottom. For some T/T c , one sees that the clusters are diffused and 
solutions are scattering to interpolate two clusters. Then, the minimum valley of x 1 seems to become flat between the 
two minima. 
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Figure 5: Temperature dependence of the residue for the vector channel peak fitted with a Breit-Wigner + continuum (left) and as Gaussian + 
continuum (right). Red points are corresponding to true minima, and gray points stand for local minima. Fit range is fixed to 7.0-12.0 GeV. 

The red point at each T/T c is the true minimum point, where x 1 hits the minimum. One sees that at low T/T c , the 
true minimum is located in the top cluster, while it moves down to the middle cluster at around T/T c ~ 1.5 — 1.6, and 
further down to the bottom cluster at around T/T c ~ 1.9 — 2.0. At T/T c > 2.5, we have only one stable solution. The 
peak positions and the continuum thresholds from the fitting with Breit-Wigner + continuum form is shown in Fig. [6] 
Both the peak positions and continuum threshold undergo a transition towards smaller values at T/T c ~ 1.5 - 2.0. 
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Figure 6: Temperature dependence of the fitting parameters for the vector channel peak fitted with a Breit-Wigner + continuum. Left : peak 
position, right : continuum threshold. Red points are corresponding to true minima, and gray points stand for local minima. Fit range is fixed to 
7.0-12.0 GeV. 



To make sure that this fitting analysis is valid, we have repeated the same procedure with an alternate fitting 
function, namely, the Gaussian + continuum form. The results are shown in the right panel of Fig. [5] The qualitative 
behaviors are the same as the previous case. It however shows only two clusters of local minima and the transition 
from the higher to the lower cluster at T/T c ~ 1.8 - 1.9. Comparing these two fitting results, we conclude that the 
qualitative behavior of the residue is independent of fitting peak function form. 

It can be concluded from the above results that the residue of T peak decreases gradually with increasing tem- 
perature and becomes a constant value at higher temperature. Especially, the rapid reduction of the residue is seen 
at T/T c — 1.5 - 2.0. It should be noted that this behavior does not directly imply that the excited states, T(2S) 
and Y(3S), disappear at lower temperature than ground state Y(1S) because one cannot eliminate other possibilities 
such as simultaneous reduction of the ground and excited states. Nevertheless, if we assume the disappearance of 
the excited states at lower temperature than the ground state, our results suggest that Y(2S ) and T(3S ) disappear at 
T/T c = 1.5 - 2.0, while Y(1S) survives up to T/T c = 3.0. 



4. Conclusion 

In summary, we have analyzed the bottomonium spectral functions at zero and finite temperature by using a newly 
developed analysis method of QCD sum rules. The Maximum Entropy Method (MEM) is adapted to extract the 
spectral function from the sum rule. 

At T — 0, the lowest peak has been obtained for each channel corresponding to Y, i]i,,Xbo, and^M . Although these 
mass spectra agree qualitatively with experimental values, their peak positions shift slightly to higher energies. By 
analyzing mock data for the vector channel and evaluating the obtained residue, we conclude that the disagreement is 
caused by the contribution of the excited states. 

Next, we have investigated the temperature dependence of the spectral functions. Temperature dependences are 
taken into account in the gluon condensates, which are estimated from the quenched lattice QCD data at finite tem- 
perature. As a result, we have found that the spectral functions of bottomonia are modified much slower as functions 
of T/T c than those of charmonia, in which the lowest peak disappears suddenly at the vicinity of T c . Using the defi- 
nitions of the upper and lower limits of the melting temperature given in Section 1331 we find that Y and r\\, survive as 
a peak in the spectral functions up to some temperature restricted to the regions of T/T c > 2.3 and T/T c > 2.1, while 
the dissociation temperatures of Xbo and^i are confined to T/T c = 1.3 - 2.5 and T/T c < 2.5, respectively. 

It should be noted, however, that our definition inevitably contains some ambiguity due to the limitation of the 
OPE and MEM. Therefore, respective results on the melting temperature should be regarded as qualitative guides. 
Furthermore, both the P-wave peaks are not found to be fully significant statistically even at T — 0, which means that 
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we cannot make a definite conclusion about the fate of these states at finite temperature. To obtain more conclusive 
results on their behavior, further studies are needed once more precise information on the OPE is available. The current 
prediction of the melting temperatures depends on the extracted temperature dependences of the gluon condensates. 
For these we have used the quenched lattice QCD data for the energy density and pressure. To go beyond the quenched 
approximation, a more detailed analysis will be required to include full QCD information on the gluon condensates, 
which will be the subject of a future investigation. 

Our results are qualitatively consistent with previous QCD sum rule analysis in conventional method 111 711 . As 
mentioned above, however, it turns out that the lowest peaks contain the excited states as well as the ground state, so 
that a deformation of such a peak depends on the behavior of the excited states. Therefore, to extract more detailed 
information on the spectral function of the vector channel, we have investigated the temperature dependence of the 
residue of the lowest peak obtained from MEM for this channel. Then, we have observed that the residue decreases 
with increasing temperature. The results are consistent with a picture that the excited states, T(2S) and Y(3S), 
dissociate at lower temperatures than the ground state T(15). Our results indicate that T(2S) and T(3S) disappear in 
the temperature region of T/T c = 1 .5 - 2.0. 
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